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异 质 生境 下 山东 银 莲花 的 转录 组 分 析 及 EST-SSR 分 子 标记 开发 


BIN, HER, JAR, KEE, KERR, FWE, TE 


(烟台 大 学 生命 科学 学 院 ， 山 东 烟台 264005) 
摘 要 : 为 探究 山东 银 莲 花 在 全 光照 的 山顶 灌 从 和 阴暗 的 针 阔 混交 林 下 两 种 不 同 生 境 中 的 生态 适应 机 种 
并 开发 其 EST-SSR 分 子 标记 ， 该 研究 利用 lumina 高 通 量 测序 技术 对 开花 期 的 山东 银 莲 花 叶 片 进行 转 


PN 


录 组 测序 ， 获 取 其 功能 注释 和 差异 表达 基因 。 结 果 表明 : 〈1) 转录 组 测序 共 得 到 53 536 条 Unigenes 序 


列 ， 其 中 27 448 条 成 功 获得 注释 。 (2) 差异 表达 基因 5 635 个 ，1 600 个 在 山顶 灌 丛 的 山东 银 莲花 中 上 


调 表 达 ， 其 余 4035 个 下 调 表 达 。 有 2 460 个 差异 表达 基因 注释 到 GO 数据 库 2 533 个 三 级 条 目 中 ，1 


051 个 差异 表达 基因 注释 到 KEGG 数据 库 的 113 条 代谢 通路 中 。 G) 山东 银 莲花 适应 于 异 质 生境 的 代 
谢 通路 主要 涉及 光合 作用 -天 线 蛋 白 通 路 和 类 黄酮 生物 合成 通路 ， 光 合作 用 -天 线 蛋 白 通 路 中 lhca5 基因 
EWR, lhcal-3 基因 下 调 表 达 ， 类 黄酮 生物 合成 通路 中 chs. c4hy Boh, Bh, flsy ans, chi, 
ccoaomt 和 het 基因 均 上 调 表 达 。 (4) 从 山东 银 莲 花 转 录 组 数据 中 共 搜索 获得 7 146 个 SSR 位 点 分 布 于 
6 006 条 Unigenes 序列 中 ， 共 计 106 种 重复 基 序 ， 优 势 重复 基 序 为 单 核 昔 酸 重复 。 设 计 合 成 100 对 EST- 
SSR 引物 中 共有 68 对 引物 具有 有 效 性 ， 其 中 11 对 具有 多 态 性 ， 共 扩 增 24 个 多 态 性 片段 。 该 研究 结果 


na 


有 助 于 更 深入 地 理解 山东 银 莲 花 在 不 同 生 境 中 的 适应 性 调节 机 制 ， 并 首次 开发 其 EST-SSR 分 子 标记 填 


补 该 方面 的 空白 ， 为 该 生物 的 保护 和 利用 提供 了 重要 的 分 子 标记 资源 。 
关键 词 ， 异 质 生境 ， 山 东 银 莲花 ， 转 录 组 ，EST-SSR， 和 生态 适应 机 币 
中 图 分 类 号 : Q943 文献 标识 码 : A 
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Transcriptome analysis and development of EST-SSR molecular markers 


in Anemone shikokiana under heterogeneous habitats 


SHAN Xiaohan, AN Kang, ZHOU Chunxia, ZHANG Xin, PANG Yujuan, LI Lixia, BIAN Fuhua 


(College of Life Sciences, Yantai University, Yantai 264005, Shandong, China ) 


Abstract: In this study, the ecological adaptation mechanisms of Anemone shikokiana in two distinct habitats, 
namely full-light hilltop scrub and shady mixed broadleaved-coniferous forest, were researched by Illumina 


high-throughput sequencing technology for leaves collected during the flowering stage. Moreover, EST-SSR 
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molecular markers were developed based on SSR locus distribution characteristics. The results were as follows: 
(1) A total of 53 536 Unigenes sequences were obtained, of which 27 448 were successfully annotated. (2) 5 
635 DEGs were obtained after filtering the low abundance genes, 1 600 up-regulated and 4 035 down-regulated 
genes comparing A. shikokiana in full-light hilltop scrub and in shady mixed coniferous forest. GO 
classification results showed that 2 460 DEGs were annotated to 2 533 tertiary entries. In addition, 1 051 DEGs 
were involved in 113 KEGG pathways. (3) The comprehensive analysis of the photosynthesis-antennal 
proteins pathway related genes revealed that the expression of /hca5 was significantly higher, while the 
expression of lhcal, lhca2 and lhca3 was significantly lower. Meanwhile, the comprehensive analysis of the 
flavonoid biosynthesis pathway related genes revealed that the expression of chs, c4h, f3’h, fh, fls, ans, chi, 
ccoaomt and hct was significantly higher. (4) A total of 6 006 unigenes containing 7 146 SSRs were identified 
among 53 536 unigenes by using MISA software from the transcriptome data of A.shikokiana. In the identified 
SSRs, the dominant repeat motifs were single nucleotide repeats in 106 repetitive motif types. Among the 100 
pairs of EST-SSR primers, a total of 68 pairs were effective and 11 pairs with polymorphism, and 24 
polymorphic fragments were amplified. Overall, in this paper, the adaptation mechanisms of A.shikokiana in 
heterogeneous habitats were analyzed at the molecular level, and for the first time EST-SSR molecular markers 
were developed to fill the gap in this area, which would provide important molecular marker resource for the 
conservation and utilization of this species. 
Key words: heterogeneous habitat, Anemone shikokiana, transcriptome, EST-SSR, ecological adaptation 
mechanisms 

山东 银 莲 花 (Anemone shikokiana) fe RF ARIE TER WN E EE AN, Ta tod i PB JG AR 
名 和 日 本 四 国 岛 〈 王 起 等 ，2014) 。 在 中 国 被 列 入 《世界 自然 保护 联盟 濒危 物种 红色 名 录 》 (IUCN) 
易 危 (VU) 等级， 为 典型 的 稀有 物种 〈 陈 春 利 等 ，2018) 。 其 根 可 入 药 ， 有 清热 解毒 、 止 血 除湿 的 功 
RM CUBR, 2014) ; 因 其 株 型 优美 、 花 大 、 花 期 长 ， 在 园艺 上 具有 潜在 的 开发 价值 ( 侯 元 同和 刘 冰 ， 
2010) 。 

山东 银 莲花 生长 在 海拔 约 600 m 以 上 的 山顶 灌 从 和 人 针 痊 混交 林 下 两 种 异 质 生 境 中 ， 山 东江 从 常年 
全 光照 ， 土 壤 是 以 碎 石 、 沙 粒 为 主 的 砂 质 性 土壤 ， 而 针 闪 混交 林 下 阴暗 潮湿 ， 土 壤 为 营养 丰富 且 朴 松 
AS (Pang et al., 2020) 。 岳 喜 元 等 2023) 研究 表明 黄河 三 角 洲 的 芦苇 通过 调整 个 体 大 小 与 叶 性 
状 适 应 异 质 生 境 下 水 分 条 件 变化 ， 王 妍 方 等 〈2023 ) 通过 木棉 表 型 可 塑性 探究 其 适应 干 热 河 谷 和 热带 
雨林 两 种 极端 环境 的 机 制 。 山 东 银 莲花 则 通过 改变 叶片 、 根 系 的 形态 结构 及 生理 等 策略 适应 完全 不 同 
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的 生长 环境 (FICK, 2019; Pang er al., 2020; EW, 2021) 。 然 而 ， 至 今 对 该 种 适应 异 质 生境 
的 分 子 机 制 尚 不 清楚 。 植 物 适 应 环境 的 本 质 是 基因 特异 性 表达 的 结果 ， 基 因 的 特异 性 表达 一 定 程度 上 
由 转录 水 平 进行 调控 〈 张 棒 雨 等 ，2007) 。 转 录 组 作为 特定 组 织 或 细胞 在 某 一 发 育 阶段 或 功能 状态 下 
所 有 mRNA 的 集合 ， 是 研究 基因 功能 与 结构 的 有 效 手 段 〈 和 鲁 艳 辉 等 ，2021) 。 运 用 Numina 高 通 量 测 
序 技术 对 植物 转录 组 进行 测序 分 析 ， 可 高 效 、 精 准 地 获得 大 量 特定 组 织 或 细胞 中 表达 的 基因 序列 ， 实 
现在 无 参考 基因 组 的 基础 上 进行 分 析 〈 郭 连 安 等 ，2021) 。 同 时 转录 组 分 析 应 用 于 两 色 金鸡 菊 〈 孙 浩 
BF, 2022). ESM mia (SAMS, 2022) 等 多 种 植物 ， 且 已 广泛 应 用 于 探究 不 同 生 境 下 植物 
的 适应 性 ( 周 春 苗 等 ，2022;， 韦 陈 彬 等 ，2022) 。 

分 子 标记 可 用 于 植物 的 系统 进化 分 析 、 功 能 基因 标记 和 分 子 辅助 育种 等 ， 随 着 高 通 量 测序 技术 的 
发 展 ， 由 转录 组 测序 得 到 的 EST-SSR 分 子 标记 被 作为 最 常用 的 分 子 标记 技术 之 一 。 其 操作 简单 ， 具 有 
< 显 性 遗传 、 稳 定性 好 和 准确 度 高 等 优 于 其 他 分 子 标记 的 特点 《后 远 濒 等 ，2022) 。EST-SSR 标记 开 
发 成 本 低 ， 适 用 于 无 参考 基因 组 的 物种 ， 可 有 效 地 推动 分 子 标记 技术 在 植物 遗传 多 样 性 分 析 、 指 纹 图 
谱 构建 、 种 质 资源 保护 、 品 种 鉴定 等 方面 的 应 用 《杨雄 等 ，2021) ， 目 前 已 在 大 麻 (边境 等 ，2023) 、 
PRL CHR AGAR, 2022) . WE FEF, 2022) 等 植物 中 开发 与 应 用 。 而 关于 山东 银 莲花 分 子 标记 
的 报道 尚 不 多 见 ， 未 见 公开 报道 的 可 利用 SSR 标记 ， 由 此 使 得 山东 银 莲 花 在 分 子 标记 方面 的 研究 难以 
深入 。 

本 研究 以 异 质 生境 下 山东 银 莲 花 为 研究 对 象 ， 利 用 Illumina 高 通 量 测序 技术 获取 其 转录 组 数据 ，3 
进行 组 装 和 功能 注释 ， 得 到 其 差异 表达 基因 和 大 量 丰 富 的 SSR 位 点 ， 随 后 进行 EST-SSR 引物 设计 并 加 
以 验证 ， 以 期 探讨 山东 银 莲 花 在 异 质 生境 中 的 生态 适应 机 制 ， 并 开发 适用 的 EST-SSR 分 子 标记 ， 填 补 
山东 银 莲 花 在 以 上 两 方面 的 研究 空白 ， 为 今后 其 遗传 多 样 性 分 析 、 种 质 资源 评价 、 功 能 基因 标记 和 分 
了 辅助 育种 等 研究 莫 定 基础 ， 进 一 步 实现 物种 保护 与 发 展 利 用 。 
1 材料 与 方法 
1.1 实验 材料 

本 研究 所 用 实验 材料 均 为 开花 期 的 山东 银 莲 花 叶 片 ， 其 中 转录 组 测序 材料 于 2021 年 6 月 取 自 山东 
省 烟台 市 昆 帘 山 国 家 级 自然 保护 区 的 泰 磷 顶 和 骞 风 岭 ( 表 1) ， 引 物 筛选 材料 于 2021 年 9 月 取 自 青岛 
崂山 和 烟台 昆 寄 山 ( 表 2) 。 采 和 集 后 用 蒸馏 水 清洗 干净 ， 立 即 放 入 液 氮 中 冷冻 。 转 录 组 测序 由 上 海 欧 易 
生物 医学 科技 有 限 公司 完成 ， 引 物 筛 选材 料 带 回 实验 室 后 采用 改良 的 CTAB 法 09 提 取 DNA， 并 利用 
1% 琼 脂 糖 凝 胶 电泳 对 所 提取 的 DNA 进行 检测 ， 将 检测 合格 的 DNA 统一 稀释 至 20 ng'hL1， 放 入 -20 °C 


冰箱 中 保存 备用 。 


表 1 用 于 转录 组 测序 的 山东 银 莲花 采样 信息 


Table 1 The Sampling information of Anemone shikokiana for transcriptome sequencing 


pa yape Tie ETEN 
ai 地 点 生境 经 度 纬度 o eee eas 
Serial : ; ; . Altitude [Illumination Temperature 
Location Habitat Longitude Latitude 3 
number (m) (dux) CC) 
HE Be at LL ae ie Tit 山顶 灌 从 
D1 Yantai Kunyu Hill b 121°78'42.7" E 37°15'26.8" N 792.0 23 600 24.4 
Mountain Tai Bo Ding atop Seru 
HE Be at LL ase Tt 山顶 灌 从 
D2 Yantai Kunyu Hill b 121°78'43.1" E 37°15'25.3" N 780.6 23 600 24.4 
Mountain Tai Bo Ding ilitop:scru 
E Be ar LL ae aT 山顶 灌 从 
D3 Yantai Kunyu Hillt b 121°78'47.5" E 37°15'25.9" N 781.4 23 600 24.4 
Mountain Tai Bo Ding Eee a 
WAG ERURKI PEPIRESEAE 下 
Yantai Kunyu axed 
L1 y 121°77'40.6" E 37°14'56.7" N 674.7 500 18.7 
Mountain Cold Wind broadleaved 
Ridge coniferous 
forest 
Yantai Kunyu 
L2 x 121°77'31.4" E 37°14'53.1" N 674.9 500 18.7 
Mountain Cold Wind proadleaved 
Ridge coniferous 
forest 
Yantai Kunyu Mixed 
L3 z 121°77'53.9" E 37°14'48.3" N 675.3 500 18.7 
Mountain Cold Wind broadleaved 
Ridge coniferous 
forest 
表 2 用 于 引物 筛选 的 山东 银 莲 花 采 样 信息 
Table 2 The Sampling information of Anemone shikokiana for EST-SSR molecular markers 
编号 地 点 经 度 纬度 海拔 Cn) 
Serial number Location Longitude Latitude Altitude 
oe N Wee MaN 
QL2 f aw Liig 120°37'50.8" E 36°11'23.6" N 962.33 
Qingdao Laoshan Slippery Pass 
HE EE, 8a LL Ze he 
YDI . AS Ee ahi i f 121°46'13.1" E 37°14'51.0" N 870.00 
Yantai Kunyu Mountain Tai Bo Ding 
bs zB a 
YLI . MA eki | 老 铁 山 121°41'49.4" E 37°16'00.1" N 582.00 
Yantai Kunyu Mountain Laotie Mountain 
T S L WF TL PL ree 
QDI Qingdao Laoshan Laodeng Scenic Area 120°37'43.3" E 36°10'39.0" N 1 019.02 
Danfeng Peak 
Æ N IE MaN 
QLI . aw Liig 120°37'50.8" E 36°11'23.6" N 962.33 
Qingdao Laoshan Slippery Pass 
1.2 转录 组 分 析 


De novo 拼接 : 使 月 


根据 序列 相似 愧 


一 套 最 终 的 Unigene， 以 此 作为 后 续 分 析 的 参考 序列 。 


功能 注释 : 利 


Prot, eggNOG, KEGG 数据 库 以 及 利 


之 后 再 利 月 


H CD-HIT 软 从 


H Trinity (Manfred G et al., 2011) 软件 paired-end 的 拼接 方法 得 到 Transcript 序列 ， 
FE 以 及 长 度 ， 挑 选 出 最 长 的 一 条 作为 Unigene， F 聚 类 去 见 余 得 到 
| diamond (Benjamin et al., 2015) 软件 将 Unigene 比 对 到 NR, KOG, GO, Swiss- 


J HMMER 软件 比 对 Pfam 数据 库 来 进行 Unigene 的 功能 分 析 。 


? 


reads 数 进行 差异 显著 性 检验 ， 以 差异 倍数 |log2FoldChanges|>1 及 差异 显 


差异 表达 基因 往 选 :利用 DESeq (Anders S et al., 2012) 软 伯 


住 化 处 理 ， 采 用 BaseMean (82% (i # 


表达 量 ， 


1.3 引物 设计 与 筛选 


越 高 则 应 有 


基于 山东 银 莲 花 转录 组 测序 结果 ， 使 
重复 基 序 的 长 度 为 1、2、3、4、5 #6, Rea 
间 的 序列 长 度 不 超过 100bp 即 被 认定 为 复合 微 卫星 。SSR 作为 一 种 用 途 广泛 的 分 子 标记 ， 


日 价值 


Jif 


F 对 各 个 样本 基 


因 的 counts 数目 进行 标 


计算 差异 倍数 ， 并 采用 NB 负 二 项 分 布 检 验 的 方式 ) 对 


著 性 检验 结果 p<0.05 为 标准 筛 


LK (ARES, 2012) 。 当 SSR 长 度 宇 20 bp 时 可 显示 出 较 高 的 多 态 拆 


F MISA (Sebastian B et al., 2017) 进行 SSR 预测 ， 设 置 


E 复 次 数 分 别 为 10、6、5、5、5 和 5 次 ， 两 个 微 卫 星之 


其 多 态 性 水 平 
E， 长 度 在 12~20 


bp 之 间 相 对 较 低 ， 在 12 bp 以 下 则 很 低 〈Temnykh et al., 2001) 。 根 据 SSR 长 度 宇 20 bp、 去 除 复 合 型 


SSR 位 点 和 容易 与 多 腺 背 化 作 月 


择 100 对 SSR 位 点 利 月 


ANS100， 送 至 北京 六 合 华 大 基因 科技 有 限 公司 合成 。 


ME. PCR 扩 增 体系 为 20 uL: 2xMaster Mix 4 uL, E] 


和 dd H2010 hL。 反 应 程序 : 94 CHAH 


个 循环 ，72 'C 终 延 作 


利用 青岛 崂山 滑 溜 口 所 采集 样品 通过 1% 琼 脂 糖 凝 胶 


可 正确 扩 增 的 引物 利用 样品 YDI ( 


L Ue TH 


方法 进行 目 


区 丹 炉 峰 附近 ) ~ QLI CE 


2 结果 与 分 析 
2.1 转录 组 数据 分 析 


2.1.1 总 体 特 征 
分 布 在 6.98~7.10 G，Q30 碱 基 分 布 帮 
续 分 析 。 利 用 De novo HHF H 
最 短 Unigene 为 301 bp， 平 均 长 度 为 1081.69 bp。 用 皮尔 逊 相关 系数 表示 D AIL A 
HI ERE) , D 组 间 生 物 学 重复 的 相关 性 均 大 了 
间 相关 性 均 大 了 


2.1.2 功能 注释 对 拼接 组 装 的 53 536 条 Unigenes 进行 常见 功能 数据 库 注 释 ， 在 7 个 数据 库 中 均 兴 


E URF 0.8 aA KM 


对 山东 银 莲花 进行 无 参 转录 组 分 析 ， 共 获得 42.34 G 的 Clean Data， 各 样本 的 有 效 数 据 量 


10 min，4'C 保 存 。 随 后 将 无 法 正 


JAG Ear Ze he) ~ YLI H 


F 0.831 8 ， 表 明 样 品 


ANS 


岛 崂 山 滑 


E 3 min; 94 ‘C4SPE 30s, 55 'C 退 火 35 s，72 'C 延 人 


Tl 


E 93.68%~94.20% 之 间 ， 证 明 本 次 测序 数据 准确 性 忆 


昌 相 混杂 的 A/T 重复 基 序 的 原则 将 所 有 SSR 位 点 进行 初步 筛选 后 随机 选 


日 软件 Primer3 (Rozen S & Skaletsky H, 2000) 进行 引物 设计 ， 命 名 为 ANSOOI- 


BE 泳 对 100 对 引物 进行 筛选 ， 以 检测 引物 有 效 


下 游 引 物 (10 umol L!) 各 2nL、DNA 模板 2 uL 
160s, 35 
Ay Tg BT PP KD E S POR, 
Gear GL) 、QD1 CFAI 
YE AL) 通过 8% 非 变性 聚 丙烯 酰胺 凝 胶 电泳 和 银 染 显 色 的 
的 条 带 的 分 离 和 检测 ， 筛 选 出 具有 多 态 性 的 EST-SSR 引物 。 


好 ， 可 用 于 后 


H Unigene 53 536 条 ， 总 长 度 为 57 909 453 bp， 最 长 Unigene 为 15 502 bp, 


昌 间 相关 愧 


fF 品 间 基 


因 表 达 水 平 


F 0.8689, 工 组 间 生 物 学 重复 之 


FE 较 高 ， 可 以 此 作为 后 续 分 析 的 参考 序列 。 


EFE 释 成 功 


的 Unigenes 有 3 505 条， 占 总 数 的 6.55%, EDE 1 个 数据 库 中 注释 成 功 的 Unigenes 有 27 448 条 ， 占 
总 数 的 51.27% (图 1) 。 在 NR 数据 库 中 成 功 注释 的 条 数 最 多 ， 共 有 27 044 条 ， 占 总 数 的 50.52%， 其 


余 分 别 为 eggNOG Æ (23 934, 44.71%) 、Swissprot JÆ (19 637, 36.68%) 、Pfam JÆ (17 859, 


33.36%) ~ GO Æ (17540, 32.76%) 、KOG JÆ (5481, 29.19%) 和 KEGG Æ (5481, 10.24%) 。 
按照 NR 数据 库 中 山东 银 莲 花 的 Unigenes 45 Fk (ty Ft DL We By 7, A ADA FE Ee es A a BH SF SE 


(Aquilegia coerulea，56.72%)， 其 次 按 匹 配 度 由 高 到 低 依次 为 博 落 回 (Macleaya cordata, 6.08%). #5} Æ 


(Papaver somniferum, 4.06%). %¥(Nelumbo nucifera，3.29%)、 和 葡萄 (Vitis vinifera，1.92%)、 欧 洲 栓 皮 
栎 (Quercus suber，0.94%)、 月 季 花 (Rosa chinensis, 0.9%), Sphaerulina musiva SO2202 (0.73%) 、 中 华 
Sik WK WK (Actinidia chinensis var. chinensis, 0.58%). iE] H 3% (Helianthus annuus，0.55%) 和 其 他 物种 


(24.24%) 。 
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上 方 条 形 图 的 数字 代表 下 方 矩 阵 中 对 应 加 黑 点 的 数据 库 交 集 的 结果 ， 左 侧 柱状 图 代表 各 数据 库 全 部 注释 到 基因 数量 。 


The numbers in the top bar represent the results of the intersection of the databases corresponding to the black dots in the matrix 


— 


below, and the bars on the left represent the number of genes annotated to each database. 


图 1 各 数据 库 注释 韦 恩 图 


Fig. 1 Venn diagram for each database annotation 
2.2 差异 表达 Unigene 分 析 
2.2.1 差异 表达 Unigene 筛选 对 山东 银 莲 花 的 转录 组 数据 进行 筛选 ， 共 挖掘 差异 表达 Unigene 5 635 个 ， 
数量 较 多 ， 其 中 与 林 下 (CL) 的 山东 银 莲花 相 比 ， 山 项 CD) 的 山东 银 莲花 中 有 上 调 基因 1 600 个 ， 下 


调 基因 4 035 个 《图 2) ， 说 明 不 同 生长 环境 对 山东 银 莲花 基因 表达 影响 较 大 。 


* 显 著 上 调 1 600 
Significant up 1 600 


* 显 著 下 调 4 035 
Significant down 4035 


不 显著 


Non-significant 


SHEP {ANY F247 HX -log,, P value 


-10 -5 0 5 10 
差异 倍数 的 对 数 log, FoldChange 


图 2 差异 表达 基因 的 火山 


Fig. 2 Volcano map of differentially expressed genes 

2.2.2 差异 表达 Unigene 富 集 分 析 将 山东 银 莲 花 的 差异 表达 Unigene 进行 GO (K3) 和 KEGG 富 集 分 析 

(图 4) ， 结 果 显 示 共 有 2 460 个 差异 表达 基因 被 注释 到 GO 数据 库 2 533 个 三 级 条 目 中 ， 共 有 1 051 个 
差异 表达 基因 被 注释 到 KEGG 数据 库 的 113 条 代谢 通路 中 。 

根据 GO 数据 库 一 级 条 目 生 物 学 过 程 (biological process, BP)、 细 胞 成 分 (cell component, CC) 和 分 子 

功能 (molecular function, MF) 的 不 同 ， 筛 选 三 种 分 类 中 对 应 差异 基因 数目 大 于 2 的 GO 二 级 条 目 ， 按 照 

每 个 条 目 对 应 的 尸 值 由 小 到 大 排序 可 得 前 10 个 富 集 分 类 ， 山 东 银 莲花 的 差异 表达 基因 主要 显著 富 集 在 


脱落 酸 代 谢 过 程 (abscisic acid metabolic process GO:0009687) 、 注 蓓 醇 生物 合成 过 程 (menthol 


biosynthetic process GO:0031525) 、 硫 脂 生物 合成 过 程 (sulfolipid biosynthetic process GO:0046506) 等 


4 


Pa 


生物 学 过 程 中 ， 胞 质 大 核糖 体 亚 基 (cytosolic large ribosomal subunit GO:0022625) 、 胞 质 小 核糖 体 亚 


(cytosolic small ribosomal subunit GO:0022627) 、 核 糖 体 (Ribosome GO:0005840) 等 细胞 成 分 中 ， 芳 


醇 脱氧 酶 (NADP+) 活性 [aryl-alcohol dehydrogenase (NADP+) activity GO:0047681]、 异 戊 烯 醇 脱 氧 酶 


活性 (isopiperitenol dehydrogenase activity GO:0018458) 、 卡 维 醇 脱氧 酶 活性 Ccarveol dehydrogenase 


activity GO:0018459) 等 分 子 功能 中 。 
将 山东 银 莲 花 的 差异 表达 基因 比 对 到 KEGG 数据 库 的 参考 代谢 通路 中 ， 可 划分 到 6 个 类 别 中 ， 其 
中 与 代谢 相关 通路 的 Unigenes 数量 最 多 ， 有 480 个 ， 占 比 45.67%; 其 次 遗传 信息 加 工 相 关 通 路 的 


下 


Unigenes 有 425 个 ， 占 比 40.44%; 与 环境 信息 处 理 相关 通路 的 Unigenes 有 58 个 ， 占 比 5.52%; 与 生物 


系统 相关 通路 的 Unigenes 有 50 个 ， 占 比 4.76%; 与 细胞 过 程 相关 通路 的 Unigenes 有 37 ^, Ek 


3.52%。 共 有 10 & KEGG 通路 达到 极 显 著 水 平 (P<0.01) ， 分 别 为 核糖 体 (ribosome ko03010)、 光 合作 


用 -天 线 和 蛋白 (photosynthesis - antenna proteins ko00196)、 类 黄酮 生物 合成 (flavonoid biosynthesis ko00941)、 


倍 半 蓝 和 三 茧 生物 合成 (sesquiterpenoid and triterpenoid biosynthesis ko00909) . 、 酮 体 的 合成 与 降解 


(synthesis and degradation of ketone bodies ko00072)、 谷 胱 甘 肽 代谢 (glutathione metabolism ko00480)、 植 


物 -病原 互 作 (plant-pathogen interaction ko04626), MAPK 信号 通路 -植物 (MAPK signaling pathway - plant 


ko04016)、 黄 酮 和 黄酮 醇 的 生物 合成 (flavone and flavonol biosynthesis ko00944)、 二 葵 乙 烯 类 、 二 芳 基 庚 


烷 类 和 姜 辣 素 的 生物 合成 (stilbenoid, diarylheptanoid and gingerol biosynthesis ko00945)。 
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图 3 差异 表达 基因 GO 富 集 Top30 气泡 图 
Fig. 3 GO enrichment of differentially expressed genes in the Top30 


活性 ，6、 硫 脂 生物 合成 过 程 ; 
亚 基 ，11、 细 胞 质 小 核糖 体 


FE 


生物 发 生 ，17、 核 糖 体 大 亚 基 组 装 ; 


18、 细 胞 质 翻 译 延伸 ; 


十 
结合 ; 


23、 


LEA 


SHAS; 


质 标 签 ， 24、 序 列 特异 性 DNA 


22、 光 系统 ; 


7、 薄 荷 醇 生 物 合成 过 程 ，8、 核 糖 体 的 结构 成 分 ，9、 细 胞 质 翻译 ; 
其 ，12、 翻 译 ; 13、 对 几 丁 质 的 反应 ; 14、 核 糖 体 |; 15、 光 合作 


酶 活性 ， 3、 胆 碱 激酶 活性 ，4、 蜡 哌 喧 醇 脱氧 酶 活性 ，5、 芳 基 醇 脱氧 酶 (NADP+) 


10、 细 胞 质 大 核糖 体 


]， 光 收集 ; 16、 核 糖 体 


19、photosystem; 20、DNA 结合 转录 


25、 细 


胞 外 区 域 ， 26、(S)- 去 


27、 柠 榜 酸 裂解 酶 复合 物 ，28、 吞 噬 泡 ;29、Pplastoglobule; 30、Set3 REW. 


1. Abscisic acid metabolic process; 2. Carveol dehydrogenase activity; 3. Choline kinase activity; 4. Isopiperitenol dehydrogenase 


天 


F 


子 活 性 ，21、 叶 绿 素 
古 柯 碱 合成 酶 活性 ; 


activity; 5. Aryl-alcohol dehydrogenase(NADP+) activity; 6. Sulfolipid biosynthetic process; 7. Menthol biosynthetic process; 8. 


Structural constituent of ribosome; 9. Cytoplasmic translation; 10. Cytosolic large ribosomal subunit; 11. Cytosolic small 


ribosomal subunit; 12. Translation; 13. Response to chitin; 14. Ribosome; 15. Photosynthesis, light harvesting; 16. Ribosome 


biogenesis; 17. Ribosomal large subunit assembly; 18. Cytoplasmic translational elongation; 19. Photosystem; 20. DNA-binding 


transcription factor activity; 21. Chlorophyll binding; 22. Photosystem; 23. Protein tag; 24. Sequence-specific DNA binding; 25. 


Extracellular region; 26. (S)-norcoclaurine synthase activity; 27. Citrate lyase complex; 28. Phagocytic vesicle; 29. Plastoglobule; 


30. Set3 complex. 
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富 集 分 数 Enrichment Score 
图 4 差异 表达 基因 KEGG 富 集 Top30 气泡 图 
Fig. 4 KEGG enrichment of differentially expressed genes in the Top30 

1、 核 糖 体 ，2、 光 合作 用 -触角 和 蛋白， 3、 黄 酮 类 化 合 物 的 生物 合成 '， 4、 倍 半 莫 和 三 莫 类 化 合 物 的 生物 合成 ; 5、 酮 体 的 
合成 和 降解 6、 谷 胱 甘 肽 代谢 ，7、 植 物 与 病原 体 的 相互 作用 ;8、MAPK 信和 号 通路 -植物 ，9、 黄 酮 和 黄酮 醇 的 生物 合成 ; 
10、 二 茶 乙 烯 类 、 二 芳 基 庚 烷 类 和 姜 酚 的 生物 合成 ，11、 维 生 素 Bo 代谢 ，12、 泛 醒 和 其 他 藻类 - 泛 醒 的 生物 合成 ，13、 
牛 磺 酸 和 低 牛 磺 酸 代谢 ，14、 共 丙 氨 酸 、 酷 氨 酸 和 色 氮 酸 的 生物 合成 15、 人 花生 四 烯 酸 代谢 ;16、 核 黄 素 代谢 ，17、 4 
物 素 代谢 ，18、 植 物 激素 信号 转 导 ; 19、 莫 类 化 合 物 骨 架 的 生物 合成 ， 20、 昼 夜 节律 -植物 ，21、 类 胡 葛 卜 素 的 生物 合成 ; 
22、 抗 坏 血 酸 和 醋酸 代谢 ，23、 丁 酸 盐 代谢 ，24、 甘 油 磷脂 代谢 ，25、 瑟 氨 酸 、 亮 氨 酸 和 异 亮 氨 酸 降解 ，26、 亚 油 酸 代 
谢 ，27、 丙 烷 、 哌 啶 和 吡啶 生物 碱 的 生物 合成 ; 28、 精 氮 酸 和 膊 所 酸 代谢 ，29、 共 丙 氨 酸 的 代谢 ，30、 不 饱和 脂肪 酸 的 
生物 合成 。 


1. Ribosome; 2. Photosynthesis-antenna proteins; 3. Flavonoid biosynthesis; 4. Sesquiterpenoid and triterpenoid biosynthesis; 5. 


tt 


Synthesis and degradation of ketonebodies; 6. Glutathionemetabolism; 7. Plant-pathogen interaction; 8. MAPK signaling pathway- 
plant; 9. Flavone and flavonol biosynthesis; 10. Stilbenoid, diarylheptanoid and gingerol biosynthesis; 11. Vitamin B6 metabolism; 
12. Ubiquinone and other terpenoid-quinone biosynthesis; 13. Taurine and hypotaurine metabolism; 14. Phenylalanine, tyrosine 
and tryptophan biosynthesis; 15. Arachidonic acid metabolism; 16. Riboflavin metabolism; 17. Biotin metabolism; 18. Plant 
hormone signal transduction; 19. Terpenoid backbone biosynthesis; 20. Circadian rhythm-plant; 21. Carotenoid biosynthesis; 22. 
Ascorbate and aldarate metabolism; 23. Butanoate metabolism; 24. Glycerophospholipid metabolism; 25. Valine,leucine and 
isoleucine degradation; 26. Linoleic acid metabolism; 27. Tropane, piperidine and pyridine alkaloid biosynthesis; 28. Arginine and 


proline metabolism; 29. Phenylalanine metabolism; 30. Biosynthesis of unsaturated fatty acids. 
2.2.3 关键 差异 表达 Unigene 信息 Ae Le h Side Hh LL Ag AR HE HE a EE RG 


3) ， 其 中 1S 个 差异 表达 Unigene 属于 光合 作用 -天 线 和 蛋白 (photosynthesis - antenna proteins) 通路 ， 分 


Ny 


别 为 1 个 上 调 基 因 Ihca5 (light-harvesting complex I chlorophyll a/b binding protein 5) ~ 7 个 下 调 基因 


Ihcb\ (light-harvesting complex II chlorophyll a/b binding protein 1) 、6 个 下 调 基 因 lhcb2 ( light- 


harvesting complex II chlorophyll a/b binding protein 2) 、1 个 下调 基因 lhcb3 (light-harvesting complex IT 


chlorophyll a/b binding protein3) ; 11 个 差异 表达 Unigene 属于 类 黄酮 生物 合成 (flavonoid biosynthesis ) 


通路 ， 分 别 为 1 个 chs (chalcone synthase) ~ 1 个 c4h (trans-cinnamate 4-monooxygenase) ~ 2 个 /3'h 


(flavonoid 3'-monooxygenase) ~ 1 个 3h (naringenin 3-dioxygenase) ~ 1 ^Ñ fls (flavonol synthase) ~ 1 


个 ans (anthocyanidin synthase) . 1 个 chi (chalcone isomerase) . 1 个 


methyltransferase) 、2 个 Act (shikimate O-hydroxycinnamoyltransferase) ， 以 上 基因 皆 为 上 调 表达 。 


表 3 山东 银 莲 花 适应 于 异 质 生境 的 差异 表达 基因 分 析 


ccoaomt (caffeoyl-CoA O- 


Table 3 Analysis of differentially expressed genes in Anemone shikokiana adapted to heterogeneous habitats 


代谢 通路 通路 基因 基因 编号 
log>FoldCh 
Metabolic pathway Pathway Gene Gene ID eal Oie aie 
K08911 lhca5 TRINITY_DN9194_ c0 gl il 1 1.032205465 


TRINITY DN24746 c3 g3 i3 2 

TRINITY DN21766 c2 g7 il 1 
TRINITY _DN24746_c3_g2 13 2 
K08912 lhcb TRINITY _DN37231_cO_gl il 2 
TRINITY _DN21766_c2_g2 i2 1 


光合 作用 -天 线 蛋 TRINITY_ DN25762 c3 gl i5 2 
Phótosynthesisantenna TRINITY_DN24746 ¢3 g12 il 2 
proteins = SS ee a 
ko00196 TRINITY DN24815 c0 g2 il 2 


TRINITY_DN22809 c2 gl il3 1 
TRINITY DN23990 c2 gl i4 2 


K08193 lhcb2 
TRINITY DN24319 c0 gl i4 2 
TRINITY DN24815 c0 gl il 2 
TRINITY DN21916 c4 g3 il 1 
K08194 lhcb3 TRINITY DN24021 c0 gl il 2 
K00660 chs TRINITY_DN16756_cO_ gl il 1 
K00487 c4h TRINITY DN17787 c0 gl il 2 
TRINITY DN23510 c2 g3 il 1 

K05280 Bh 
TRINITY DN23510 ¢2 g4 il0 1 
R K00475 h TRINITY DN9341 c0 gl il 2 

类 黄酮 生物 合成 B g 

Flavonoid biosynthesis K05278 flis TRINITY _DN5340 c0 gl il 2 
ko00941 K05277 ans TRINITY DN30744 c0 gl il 1 
K01859 chi TRINITY DN25718 c0 gl il 1 
K00588 ccoaomt TRINITY DN20298 c0 gl il 1 
TRINITY DN23267 c0 g2 i2 1 

K13065 het 
TRINITY _DN21518 c0 gl il 1 

2.3 SSR 位 点 分 析 与 引物 开发 


-1.12845721 
-1.721035581 
-1.952942905 
-3.843048888 
-4.048259712 
-5.656930902 

-5.51632584 
-1.900257362 
-1.941378951 
-2.243240333 
-2.732165837 
-1.242203266 
-1.545454537 
-1.237597641 

5.944109472 

2.458589689 

8.323553071 

1.523967569 

2.813458873 

2.669799151 

2.54154359 

4.311335927 

2.677118483 

2.292870935 

1.070509174 


2.3.1 SSR 位 点 总 体 特点 对 山东 银 莲 花 转 录 组 的 53 536 条 Unigenes 序列 进行 搜索 ， 


kL 


Unigenes 序列 中 含有 7 146 个 SSR 位 点 ， 


LE HH 6 006 


条 


F 均 分 布 距离 为 8.1 kb，SSR 长 度 在 10~225 bp 之 间 ，SSR 出 


现 频率 (=SSR 个 数 / 总 Unigene 数 ) 为 13.35%，SSR 发 生 频 率 (=£ SSR 的 Unigene 数 /总 Unigene 数 ) 


为 11.22%。 山 东 银 莲花 含有 2 个 及 以 上 的 SSR 位 点 的 序列 有 907 条 ， 占 全 部 序列 的 1.68%, BARA 


重复 类 型 的 SSR 位 点 的 序列 有 420 条 ， 占 全 部 序列 的 0.78%。 山 东 银 莲花 SSR 类 型 较为 丰富 ， 共 计 


106 种 重复 基 序 


表 4) ， 优 势 重复 基 序 为 单 核 苷 酸 重复 (3248, 45.59%) ， 其 次 分 别 为 三 核 苷 酸 重复 


(2041, 28.56%) . ZERE (1 645, 23.02%) 、 六 核 苷 酸 重复 (101, 141%) 、 四 核 苷 酸 重复 


(81, 1.13%) 


. RARER (20, 0.28%) (图 $) 。 


表 4 山东 银 莲 花 中 不 同 SSR 重复 基 序 的 出 现 频率 


Table 4 Frequency of different SSR repeat motifs in Anemone shikokiana 


SSR 类 型 重复 基 序 种 类 数 重复 基 序 类 型 数量 发 生 频 率 

SSR type Number of repeat motif types Repeat motif type Number Frequency of occurrence (%) 

AIZE 7 A/T 3208 44.89 

Mono-nucleotide C/G 50 0.70 

AG/CT 927 12.97 

ZEER j AC/GT 385 5.39 

Di-nucleotide AT/AT 328 4.59 

CG/CG 5 0.07 

AAG/CTT 648 9.07 

ATC/ATG 338 4.73 

ACC/GGT 307 4.30 

AAC/GTT 205 2.87 

ZBER i AGC/CTG 188 2.63 

Tri-nucleotide AAT/ATT 168 2.35 

AGG/CCT 91 1.27 

ACT/AGT 39 0.55 

CCG/CGG 33 0.46 

ACG/CGT 24 0.34 

AAAT/ATTT 37 0.52 

IZER j7 AAAC/GTTT 10 0.14 

Tetra-nucleotide ACAT/ATGT 7 0.10 

其 他 Others 27 0.38 

AGAGG/CCTCT 4 0.06 

MEE i AAAAC/GTTTT 3 0.04 

Penta-nucleotide AAAAT/ATTTT 3 0.04 

其 他 Others 10 0.14 

ACCAGC/CTGGTG 5 0.07 

NIZE g AGATGG/ATCTCC 4 0.06 

Hexa-nucleotide AAAATC/ATTTTG 3 0.04 


其 他 Others 89 1.25 


重复 数量 Number of repetitions 


横 坐 标 为 SSR 重复 单元 长 度 ， 纵 坐标 为 该 重复 单元 的 重复 次 数 下 SSR 数量 。 


4000 


3000 
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l 2 3 4 ə 


SSR 类 型 SSR type (bp) 


The horizontal coordinate is the length of the SSR repeat unit and the vertical coordinate is the number of SSRs under the number 


of repetitions of that repeat unit. 


2.3.2 EST-SSR 引物 
EST-SSR 引物 进行 检验 ， 


图 5 山东 银 莲花 SSR 的 类 型 分 布 


Fig. 5 SSR type distribution of Anemone shikokiana 


性 聚 丙烯 酰胺 凝 胶 1 
崂 项 景区 丹 炉 峰 附 近 (QDD . FA 


引物 均 


法 ， 可 


能 在 上 述 材料 


通过 简单 的 PCR 扩 增 效果 验证 得 出 ESTSSR 引物 的 稳定 性 、 可 重复 怡 


的 开发 与 筛选 利用 1% 琼 脂 糖 凝 胶 电泳 对 从 山东 银 莲花 转录 组 数据 中 挑选 的 100 对 


可 得 68 对 有 效 扩 增 引物 。 为 探究 68 对 EST-SSR 引物 的 多 态 性 ， 采 
电泳 检验 来 源 于 烟台 昆 帆 山 泰 确 顶 (YD1) 、 烟 台 昆 宜山 老 铁 山 (YL1) 、 
BLY CQL) 的 山东 银 莲 花 DNA 扩 增 结果 ，68 对 EST-SSR 
| 中 扩 增 得 到 清晰 条 带 ， 其 中 有 11 对 EST-SSR 引物 具有 多 态 性 ( 表 5) ， 
个 多 态 性 片段 。 以 上 结果 表明 ， 基 于 山东 银 莲花 转录 组 数据 开发 EST-SSR 引物 是 一 


K 5 山东 银 莲 花 11 对 EST-SSR 引物 信息 及 扩 增 结果 


Table $ Information and amplification results of EST-SSR primers of 11 pairs of Anemone shikokiana 


] 8% 非 变 


青岛 崂山 


k 扩 增 24 


经 济 、 可 行 的 方 
E 和 多 态 性 。 


Pik Pik = 
ak 多 态 性 条 带 
名 称 及 带 数 数 
; FETC 引物 序列 Length Number 
Primer ; s Number of 
Motif Primer sequence of 1 hi 
name product amplifie po = i 
(bp)  dbands eau 
ai (A)57  CCAGTTTGGCTATGGGCTGA/GAAGGTGACGGAGCAATGGA 233 4 3 
pa (TGA)7 CCGGTTTCAGAAGCAAGCAC/CCCCTCCGGCAATTCTTTCT 239 4 2 
R (D20 TTCCCTCCTTGAACGTCTGC/AGGCAGTGGAAAAGCTCTCC 245 7 3 
oe (GA)IS | GITGATCCGTTCGTACCAGCA/CTTTGGGAAGAACGTTGGCG 122 3 1 
ANS05 (TTO)10 + TCTAGCCTGAGACGGAAGCT/ACAGACTCCAAAGTAGGCGC 125 4 2 


9 


eas (A)21 CAGTTGTCTGTGCGTGCAAA/GCTTCACGCGTCGGTTTATC 172 4 2 
ca (GA)21 GAAACCCGAAGAAGGCCAGA/ATGGATGGCTTGACGGACTG 181 4 2 
ANSO7 CAGCAGCAGGAGGAAGACAA/GACGGCGACCATCAAACAA 

n (GGA)7 A 193 4 2 
SE (GA)10 CAGCCATGAAGCATGCCATC/ATTGAAACGCCCGACTTCCT 211 3 1 
ANSO08 GCAGAGGAGGAGGAGGAAGA/CCACGACCCCCAAACTGAT 

: (TGG)7 S 226 2 1 
oe (CT)14 CGCCTAATGCTAAGCCGGTA/AGGCAGACTACAACAACGGG 163 5 5 
3 讨论 与 结论 
3.1 异 质 生 境 中 山东 银 莲 花 的 适应 机 制 分 析 

高 通 量 测序 技术 能 全 面 快速 地 获得 植物 的 转录 本 序列 信息 ， 为 植物 基因 发 据 和 次 生 代 谢 物 调 控 等 
研究 提供 了 巨大 便利 (Lu etal., 2010) 。 本 研究 利用 高 通 量 测序 技术 对 异 质 生境 中 的 山东 银 莲 花 进行 转 
录 组 测序 ， 通 过 对 其 差异 表达 基因 进行 富 集 分 析 ， 发 现 山 东 银 莲花 适应 两 种 异 质 生境 的 机 制 主要 表现 
在 光合 作用 和 类 黄酮 生物 合成 两 方面 。 

光合 作用 作为 植物 生长 和 一 切 代谢 活动 的 生理 基础 ， 受 光照 、 水 分 、 温 度 等 不 同 环境 因子 的 影响 


E 其 调控 过 程 


(Fang et al., 2022) , # 


Ht AW I MARS a/b 


chlorophyll a/b binding protein, LHC) 发 挥 重要 作 月 


的 作 月 


MS ob 
¥ He 
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基因 


大 且 可 同时 


WARE 


Zz 


含量 最 丰富 的 


faa 


Ho LHC 由 核 基 因 
日 ， 可 以 促进 光合 作 月 
家 族 主要 包含 lhca 和 lheb 两 个 进化 类 群 ， 前 者 编码 植物 光 系 统 工 中 
者 获 光 谱 中 不 同 波长 的 光 〈Mozzo et al., 2010) ， 后 
体 膜 上 50% 的 色素 形成 叶绿体 类 襄 体 


者 编 


Hk A KE 
结合 和 蛋 


(lhe) 编 


码 ， 


PERRE 


(EZE, 


本 研究 中 发 现 ， 相 比 针 阔 混交 林 下 生长 的 山 


FLAX 


Ganeteg 等 (2004) fi BA head 在 强 


阻止 叶绿体 接收 过 多 光 能 ， 


因 表 达 量 的 减少 导致 番茄 叶绿素 含量 降低 以 响应 3 
蛋白 的 合成 ， 减 弱 光 合作 用 ， 减 少 水 分 的 消耗 ， 从 而 降低 干旱 对 
丛 与 针 冰 混交 林 下 比较 ， 光 照 强 、 土 壤 湿度 小 是 最 明显 的 不 利 了 
的 上 调和 下 调 很 巧妙 的 解决 了 这 


研究 表明 ，1ca 和 Inch 两 类 基因 


实现 非 生物 胁迫 响应 。 


光 条 件 下 表达 水 3 


山顶 的 山东 银 莲 花 中 WUcb1、Whcb2 和 1hcb3 显著 下 调 ， 蒋 


东 银 莲花 ， 山 顶 灌 从 


的 山东 和 


2 & 


Á (light-harvesting complex II 
具有 迅速 将 光 能 传 到 
进行 〈Helena et al., 2023) . lhc 
的 LHC 1I 膜 蛋白 ， 捕 光 截 面 
码 植物 光 系统 本 中 的 LHC IEE 


较 


2020) 。 在 


ELE hcas 显著 上 调 ， 与 


著 提高 的 研究 结果 一 致 。 同 时 本 看 


E (2019) 研究 发 现在 干旱 胁迫 下 水 并 的 lhcb 


使 光合 作用 减弱 ，Rahele 等 2022) 的 研究 也 


FEM, 


AN 


HH lhcb 基 


F 旱 胁迫 ， 因 此 可 推测 Wcb1-3 的 下 调 可 减少 捕 光 天 线 
莲花 自身 产生 的 影响 。 山 顶 灌 


Ly AR 4 
山东 银 莲 花生 存 的 逆境 因 
难题 ， 适 应 了 山顶 灌 从 不 利 的 环境 


类 黄酮 是 植物 体内 广泛 存在 的 次 生 代 谢 物质 ， 可 吸收 280~315nm 波长 的 紫外 光 ， 保 护 植物 体 器 官 ， 


尤其 是 光合 组 织 免 受 或 少 受 辐射 伤害 〈 邹 凤 莲 等 ，2004) ， 在 植物 生长 发 育 和 热 胁 迫 等 非 生物 胁迫 响 


应 过 程 中 发 挥 多 种 生理 功能 (Chen et al., 2023; Guo et al., 2023; Zhuang et al., 2023) 。 类 黄酮 生物 合 
成 的 前 体 物 质 来 自 于 葵 丙 烷 代 谢 ， 并 受 一 系列 关键 酶 的 调控 ， 如 查 尔 酮 合成 酶 CCHS) 、 查 尔 酮 异 构 酶 
(CCHI) 、 黄 酮 醇 合 成 酶 (FLS) 、 黄 酮 3- 羟 基 化 酶 (F3H) 、 黄 酮 3' -HUE (F3'H) 、 花 色素 合成 


a 


酶 (CANS) 等 (Shen etal., 2022; RRHAYH, 2022; Xiong etal., 2016) 。 本 研究 发 现 山顶 灌 从 的 山东 银 


莲花 中 编码 相关 酶 的 基因 均 较 林 下 显著 上 调 ， 当 韵 铮 等 〈2020) 研究 表明 chs 和 chi 基因 表达 水 平 与 类 
黄酮 含量 呈 显 著 正 相 关 关系 ， 表 明山 顶 灌 从 中 的 山东 银 莲 花 在 生长 过 程 中 可 以 合成 较 多 的 类 黄酮 。 山 
顶 治 从 中 的 强 光 照 ， 除 了 可 见 光 ， 还 包括 紫外 光 ， 同 时 伴随 着 午间 的 高 温 ， 长 时 间 的 紫外 光照 射 会 造 
成 植物 DNA 及 其 它 结构 的 破坏 (Batschauer, 1993) ， 加 上 高 温 的 影响 ， 对 植物 将 造成 重大 伤害 甚至 死 
亡 。 本 研究 山顶 灌 从 山东 银 莲花 大 量 类 黄酮 的 合成 ， 可 吸收 多 余 的 紫外 光 ， 减 少 紫 外 光 对 植株 的 破坏 ， 
同时 缓解 短 时间 的 高 温 ， 以 适应 不 良 的 生活 环境 。 
3.2 山东 银 莲花 的 SSR 位 点 分 析 与 引物 开发 

山东 银 莲 花 的 SSR 类 型 相对 丰富 ， 从 单 核 苷 酸 重 复 序 列 到 六 核 苷 酸 重复 序列 不 等 。 研 究 表明 三 核 


芽 酸 基 序 中 AAG/CTT 是 主要 的 基 序 (Xiang et al., 2023) ， 与 本 文山 东 银 莲花 SSR 位 点 分 析 结 果 一 致 ， 


说 明 ESTSSR 在 其 发 生 和 进化 过 程 中 具有 高 度 保守 性 。EST-SSR 分 子 标 记 源 于 功能 基因 表达 序列 ， 具 


有 保守 程度 高 、 种 间 及 属 间 通 用 性 好 等 优点 〈 姜 永峰 等 ，2023 ) 。 本 研究 基于 拼接 质量 较 高 的 转录 组 
数据 ， 利 用 SSR 位 点 设计 并 合成 100 对 引物 ， 其 中 有 68 对 引物 可 以 实现 有 效 扩 增 ， 有 11 对 引物 在 不 
同 地 域 的 山东 银 莲花 材料 中 表现 出 多 态 性 ， 具 有 多 态 性 的 EST-SSR 分 子 标记 可 用 于 遗传 多 样 性 分 析 与 
分 子 标 记 辅 助 育种 。 因 国内 山东 银 莲 花 材 料 稀缺 、 获 取 日 本 四 国 岛 的 山东 银 莲 花 材 料 受 限 ， 在 本 研究 
中 未 进行 相关 聚 类 分 析 以 及 系统 进化 树 的 构建 ， 后 续 将 通过 银 莲 花 属 植物 对 本 研究 所 得 的 山东 银 莲 花 
EST-SSR 多 态 性 引物 的 通用 性 进行 验证 ， 同 时 利用 挖掘 的 差异 表达 基因 设计 EST-SSR 引物 用 于 区 分 不 
同 生境 的 山东 银 莲 花 居 群 ， 为 其 进化 研究 和 种 质 资源 保护 黄 定 基础 。 

综 上 所 述 ， 本 研究 对 山顶 灌 丛 和 针 疗 混交 林 下 山东 银 莲 花 的 转录 组 进行 分 析 ， 表 明 其 通过 调控 光 
合作 用 途径 和 类 黄酮 生物 合成 途径 涉及 的 基因 以 适应 异 质 生境 ， 并 首次 开发 多 态 性 较 好 的 EST-SSR 分 
子 标记 ， 填 补 山 东 银 莲花 在 此 方面 的 空白 ， 对 其 保护 与 利用 具有 重要 意义 。 
参考 文献 : 

WH, Ek, Bt, SE, 2023. KPR EST-SSR 遗传 结构 分 析 及 指纹 图 谱 构 建站 .植物 遗传 资源 学 报 : 


1 - 16[2023-10-20].https://doi.org/10.13430/j.cnki.jpgr.20230531001. [BIAN J, WANG XN, CAO K, et al., 
2023. Genetic structure analysis and fingerprint construction of hemp by EST-SSR analysis[J]. J Plant Genet 


Resour: 1—16[2023-10-20]. https://doi.org/10.13430/j.cnki.jpgr.20230531001.] 


eel, ME, XIE, AE, 2018. 山东 银 莲 花 组 织 培养 植株 再 生体 系 的 建立 中. 江苏 农业 科学 ，46(24): 


47-50. [CHEN LC, GU DZ, LIU Y, et al., 2018. Establishment of tissue culture regeneration system of 
Anemone shikokiana[J]. Jiangsu Agric Sci, 46(24): 47-50. ] 


BHA, 2022. 桑树 类 黄酮 合成 途径 关键 基因 CHIL M F3H 的 分 子 克 隆 及 功能 分 析 [D]. 镇 江 : 江苏 科技 


KÆ. [DAI MJ, 2022. Molecular cloning and functional analysis of CHIL and F3H in flavonoid biosynthesis 


pathway in mulberry[D]. Zhenjiang: Jiangsu Univ of Sci and Technology. ] 
FEZ, Kiki, WI, $, 2021. 不 同 产 地 白术 转录 组 差异 表达 基因 分 析 [J]. 中 药材 ，44(12): 2787- 


2792. [GUO LA, MO RY, TAN J, et al., 2021. Analysis on differential expressed genes in the transcriptome 


of Atractylodes macrocephalafrom different producing areas[J]. Chin Med Mat, 44(12): 2787-2792. | 


e 7c), XK, 2010. Wh A 4R 4 J. 百科 知识 ，(12) : 40. [HOU YT, LIU B, 2010. Anemone 


shikokiana(Makino) Makino[J]. Encyclopedic Knowl, (12): 40.] 


H 2019. 7K ZEA ZA EE EAE IP A AIO TD]. 南京 : 南京 农业 大 学 . [JIANG Q, 2019. 


Genomic analysis and response factors under abiotic stress in Oenanthe javanica[D]. Nanjing: Nanjing 
Agricultual University. ] 


Ai, JOR, REAR, F, 2023. 基于 转录 


ANS 
I 


测序 的 枫 香 EST-SSR 引物 开发 及 有 效 性 评价 四. 广西 植 


物 ，43(2): 327-335. [LI H, FENG YH, TANG SS, et al., 2023. Development and validity evaluation of 


Liquidambar formosana EST-SSR primers based on transcriptome sequencing[J]. Guihaia, 43(2): 327-335. ] 


XJR, 2014. 山东 银 莲 花 的 遗传 多 样 性 研究 [D]. 济南 :山东 师范 大 学 . [LIU Q, 2014. Research on genetic 
diversity of Anemone shikokiana[D]. Jinan: Shandong Normal University.] 
Zeki, HR, Amida, SE, 2022. 黄 秆 乌 哺 鸡 竹 转录 组 EST-SSR 分 子 标记 开发 与 应 用 加. 植物 科学 学 


报 ，40(3): 355-364. [LOU YF, GAN R, ZHU CL, et al., 2022. Development and application of EST-SSR 


markers based on transcriptome of Phyllostachys vivax McClure f. aureocaulis N. X. Ma[J]. Plant Sci J, 
40(3): 355-364. ] 
SHH, E RR, Wa, F, 2021. AA ATC DARA ae RRAN]. 植物 保护 学 


报 48(6): 1303-1309. [LU YH, YU FB, GUO JH, et al., 2021. Analysis of transcriptomic differences 


between winged and winglessbird cherry-oat aphid Rhopalosiphum padi on wheat[J]. Acta Phytophy Sin, 
48(6): 1303-1309.] 
逢 玉 娟 ， 孙 维 ， 王 行 松 ， 等 ，2022. 土壤 异 质 性 对 山东 银 莲 花 根 结构 的 影响 [J]. 烟台 大 学 学 报 ，35(1): 


28-34. [PANG YJ, SUN W, WANG XS, et al., 2022. Influence of soil heterogeneity on root structure of 


Anemone shikokiana[J]. J Yantai Univ, 35(1): 28—34.] 


孙 浩 男 ， 刘 冬 云 ， 李 明 阳 ， 等 ，2022. SEP mE EMF PY EB eS ee RH TI]. 分 子 植物 育种 : 


1 - 18[2023-10-20].http://kns.cnki.net/kcms/detail/46.1068.S.20221017.1205.010.html. [SUN HN, LIU DY, 
LI MY, et al., 2022. Transcriptome data analysis of Coreopsis tinctoria based on high-throughput 
sequencing[J]. Mol Plant Breed: 1—18[2023-10-20]. 
http://kns.cnki.net/kcms/detail/46.1068.S.20221017.1205.010.html.] 


HPR, Bee, FY, SB, 2022. 水 涝 胁迫 下 辣椒 转录 组 特征 分 析 及 EST-SSR tric FF K[I/OL]. 广西 植 


物 : 1 - 14[2023-10-20].http://kns.cnki.net/kems/detail/45.1134.q.20221227.1548.004.html.[TIAN HZ, 
GUO H, TIAN H, et al., 2022. Transcriptome characteristics analysis and EST-SSR marker development of 
capsicum annuum under waterlogging stress[J/OL]. Guihaia: 1—14.[2023-10- 
20].http://kns.cnki.net/kcms/detail/45.1134.q.20221227.1548.004.html. ] 


王 妍 方 ， 朱 润 军 ， 程 希 平 ，2023. 异 质 生 境 下 不 同 径 级 木棉 的 表 型 可 塑性 研究 [四 . 生态 科学 ，42(3): 61- 


66. [WANG YF, ZHU RJ, CHENG XP, 2023.Research on phenotypic plasticity of different diameter class 


Bombax ceiba in heterogeneous habitats[J]. Ecol Sci, 42(3): 61-66. ] 


EB ES, 2020. = RK Lhe 基因 的 鉴定 和 表达 模式 研究 [D]. 武汉 : 华中 农业 大 学 . [WANG YH, 2020. 


Identification and expression profiling of Lhc genes in Ramie[D]. Wuhan: Huazhong Agricultural University. ] 


EB, ERG, XIE, SE, 2014. 稀有 植物 山东 银 莲 花 (4nemone shikokiana (Makino) Makino) 的 分 布 格 


Fay Ke HE Me) A] A) HT]. 植物 研究 344): 440-445. [WANG Z, PANG YJ, LIU CL, et al., 2014. 


Distribution pattern and analysis of influencing factors of rare plant Anemone shikokiana (Makino) 


Makino[J]. Bot Stud, 34(4): 440-445. ] 


BMW, OBL, FRR, SE, 2022. SET a UU AA EA EAR ES RE AA TS]. 中 国 现 


代 中 药 ，24(12): 2342-2349. [WEI CB, GU W, TIAN R, et al., 2022. Transcriptomic data analysis of Alisma 
orientale tubers grew in different ecological environments by high-throughput sequencing[J]. Mod Chin Med, 
24(12): 2342-2349. ] 


GAGS, PAA, ASH, 2022. ET ER AE AY EE a SSR FPP TE OP PTT]. 西北 农林 科技 大 


学 学 报 〈 上 自然 科学 版 ) 50(5): 1 - 8. [WU RN, XU B, SHI FL, 2022. Analysis of SSR characteristics for 


Medicago ruthenica “Zhilixing” based on transcriptome data[J]. J Northwest A & F Univ (Nat Sci Ed), 
50(5): 1 - 8.] 
Pre, BEIE, Ri, $, 2020. 滇 黄精 转录 组 测序 及 类 黄酮 合成 相关 基因 的 分 析 [J]. 河南 农业 大 学 


学 报 ，54(6): 931-940. [XIAO YZ, HAN SM, QIN Z, et al., 2020. Analysis of transcriptome sequencing 


and related genes of flavonoids biosynthesis from Polygonatum kingianum[J]. J Henan Agric Univ, 54(6): 


931-940.] 


闫 苗 苗 ， 魏 光 成 ， 潘 效 红 ， 等 ，2008. 一 种 适用 于 动物 与 植物 总 DNA 提取 的 方法 一 一 改良 CTAB 法 中. 


安徽 农业 科学 ，(20): 8488+8558. [YAN MM, WEI GC, PAN XH, et al., 2008. A Method Suitable for 


Extracting Genomic DNA from Animal and Plant Modifed CTAB Method Agricultural[J]. J Anhui Agric 
Sci, (20): 8488+8558.] 


tat, tay, ATE, &, 2021. 白 皮 松 EST-SSR 分 子 标记 的 开发 及 应 用 四. 北京 林业 大 学 学 报 ，43(7): 


111. [YANG X, YANG N, YUAN QH, et al., 2021. Development and application of EST-SSR molecular 
markers in Pinus bungeana[J]. J Beijing For Univ, 43(7): 111.] 
于 文英 ， 高 燕 ， 逢 玉 娟 ， 等 ，2019. 山东 银 莲 花 叶 片 形态 结构 对 异 质 生境 和 海拔 变化 的 响应 四 . 生态 学 ， 


39(12): 4413-4420. [YU WY, GAO Y, PANG YJ, et al., 2019. Response of leaf morphology and structure of 
Anemone shikokiana to heterogeneous habitats and altitudechanges[J]. Acta Ecol Sin, 39 (12): 4413-4420. ] 


后 喜 元 ， 岳 伟 ， 严 宽 ，2023. BEIM = fA By Be ER 1 IRD 2k op ER PE]. 东北 林业 大 学 学 报 ， 


51(3): 18-21. [YUE XY, YUE W, YAN K, 2023. Adaptation of functional traits of Phragmites australis to 
different water Habitats in the yellow river delta[J]. J NE For Univ, 51(3): 18—21.] 


fei, MRA, WE, A, 2022. 基于 全 长 转录 组 测序 的 南京 机 EST-SSR 标记 开发 [中 . 分 子 植物 育种 : 


1-8[2023-10-20]-http://kns.cnki.net/kcms/detail/46.1068.S.20220715.1525.008.html. [YUE YH, YAN LJ, 
HUANG X, et al., 2022. Development EST-SSR markers based on transcriptome sequences of Tilia 
miqueliana Maxim[J]. Mol Plant Breed: 1—18[2023-10- 
20].http://kns.cnki.net/kcms/detail/46.1068.S.20220715.1525.008.html.] 

IKEN, Jet, GE, SF, 2007. 植物 基因 在 转录 水 平 上 的 调控 及 其 生物 学 意义 [可 . 遗传 ，(7): 793-799. 


[ZHANG CY, LONG Y FENG J, et al., 2007. Transcriptional regulation of plant genes and its significance in 


biology[J]. Hereditas, (7): 793-799. ] 


PH, BAH, RAA, 2012. 微 卫星 序列 完整 性 与 多 态 性 的 关系 四 . 基因 组 学 与 应 用 生物 学 ，31(6): 


587-591. [ZHENG Y, JIN GL, WU WR, 2012. Relationship between sequence completeness and 
polymorphism of microsatellites[J]. Genom Appl Biol, 31(6): 587—591.] 


周 春 苗 ， 蒋 莉 萍 ， 董 刚强 ， 等 ，2022. 基于 转录 组 测序 研究 不 同 海拔 生境 神农 香 菊 的 基因 表达 差异 []. 


世界 科学 技术 -中 医药 现代 化 ，24(4): 1356-1368. [ZHOU CM, JIANG LP, DONG GQ, et al., 2022. 
Analysis of differential expression genes of Chrysanthemum indicum var. aromaticum associated with 
different altitude habitats using RNA-seq technology[J]. World Sci Technol-Mod Chin Med, 24(4): 1356- 


1368.] 


BAE, FRK, MO, Se, 2004. 类 黄酮 化 合 物 在 植物 胁迫 反应 中 作用 的 研究 进展 加 . 中 国 细 胞 生 


物 学 学 报 ，(1): 39-44. [ZOU FL, SHOU SY, Ye ZZ, et al., 2004. Advances in the research on flavonoid 


biosynthesis and plant stress response[J]. Chin J Cell Biol, (1): 39-44.] 

BATSCHAUER A, 1993. A plant gene for photolyase: an enzyme catalyzing the repair of UV-light-induced 
DNA damage[J]. Plant J, 4(4): 705-709. 

BUCHFINK Benjamin, CHAO Xie, DANIEL H Huson, 2015. Fast and sensitive protein alignment using 
DIAMOND[J]. Nat Methods, 12: 59-60. 

CHEN YY, LU HQ, JIANG KX, et al., 2023. The Flavonoid Biosynthesis and Regulation in Brassica napus: A 
Review[J]. Int J Mol Sci, 24, 357. 

FANG L, MARTRE P, JIN K, et al., 2022. Neglecting acclimation of photosynthesis under drought can cause 
significant errors in predicting leaf photosynthesis in wheat[J]. Glob Change Biol, 29, 505-521. 

GANETEG U, KIMMEK F, JANSSON S, 2004. Lhca5-an LHC-type protein associated with photosystem I[J]. 
Plant Mol Biol, 54(5): 641-51 

GHANBARI MSR, TOHIDFAR M, AHMADIKHAH A, 2022. Gene profile analysis and molecular- 
physiological evaluation of tomato genotypes under drought stress[J]. J Appl Biotechnol Rep, 9(3): 747-62. 

GUO LA, TAN J, DENG XS, et al., 2023. Integrated analysis of metabolome and transcriptome reveals key 
candidate genes involved in flavonoid biosynthesis in Pinellia ternata under heat stress[J]. J Plant Res, 
136(3): 359-369. 

HELENA S, MAKIO Y, ATSUSHI T, et al., 2023. Reversible down-regulation of photosystems I and II leads to 
fast photosynthesis recovery after long-term drought in Jatropha curcas[J]. J Exp Bot, 336-351. 

LUCINSK R, JACKOWSKI G, 2006. The structure, functions and degradation of pigment-binding proteins of 
Photosystem II[J]. Acta Biochim Pol, 53 (4): 693-708. 

MANFRED G, BRIAN J, MORAN Y, et al., 2011. Trinity: reconstructing a full-length transcriptome without a 
genome from RNA-Seq data[J]. Nat Biotechnol, 29(7): 644. 

MOZZO M, MANTELLI M, PASSARINI F, et al., 2010. Functional analysis of Photosystem I light-harvesting 
complexes (Lhca) gene products of Chlamydomonas reinhardtii[J]. Biochim Biophys Acta, 1797: 212-221. 
PANG YJ, LI LX, BIAN FH, 2020. Photosynthetic and ultrastructural adaptability of Anemone shikokiana 

leaves to heterogeneous habitats[J]. Braz J Bot, 43, 979-988. 
Rozen S, Skaletsky H, 2000. Primer3 on the WWW for general users and for biologist programmers[J]. 


Methods Mol Biol, 132: 365-86. 


SEBASTIAN B, THOMAS T, THOMAS M, et al., 2017. MISA-web: a web server for microsatellite 
prediction[J]. Bioinform, 33(16): 2583—2585. 

SHEN N, WANG TF, GAN Q, et al., 2022. Plant flavonoids: Classification, distribution, biosynthesis, and 
antioxidant activity[J]. Food Che, 132531. https://pubmed.ncbi.nlm.nih.gow/35413752/. 

SIMON A, WOLFGANG H, 2012. Differential expression of RNA-Seq data at the gene level [J]. 
EMBL. https://www.bioconductor.org/help/course-materials/20 12/BioC2012/DESeq.pdf. 

STORF S, STAUBER EJ, HIPPLER M, et al., 2004. Proteomic analysis of the photosystem I light-harvesting 
antennain tomato (Lycopersicon esculentum)[J]|. Biochemistry-USA, 43(28): 9214—9224. 

TEMNYKH S, DECLERCK G, LUKASHOVA A, 2001. Computational and experimental analysis of 
microsatellites in rice (Oryza sativa L.): Frequency, length variation, transposon associations, and genetic 
marker potential[J]. Genom Res, 11(8): 1441-1452. 

XIANG NY, LU BJ, YUAN T, et al., 2023. De Novo transcriptome assembly and EST-SSR marker 
development and application in Chrysosplenium macrophyllum[J]. Genesis, 14: 279. 

XIONG S, TIAN N, LONG JH, et al., 2016. Molecular cloning and characterization of a flavanone 3- 
hydroxylase gene from Artemisia annua L.[J]. Plant Physiol Biochem, 105(5): 29-36. 

ZHUANG WB, LI YH, SHU XC, et al., 2023. The classification, molecular structure and biological 


biosynthesis of flavonoids, and their roles in biotic and abiotic stresses[J]. Molecules, 28(8): 3599. 


